Numerical simulations of reduced fluid models describing magnetized plasmas have previously been used to study the boundary of fusion plasmas extensively. A common approximation applied to these models is the so-called Boussinesq approximation. One of the requirements for the Boussinesq approximation to be applicable is that the background particle density is significantly higher than the fluctuations ...